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1.  INTRODUCTION  AND  SUMMARY 

In  this  paper  we  attempt  to  give  a  constructive,  affirma¬ 
tive  answer  to  each  of  the  following  questions. 

1.  Given  a  function  f  and  an  interval  I.  is  it  poss¬ 
ible  to  tell  a  priori  whether  or  not  one  can  accurately 
approximate  f  via  a  low  degree  rational  function? 

2.  Can  such  a  rational  function  be  easily  constructed 
explicitly,  so  that  one  encounters  no  poles  on  the  interval 
of  approximation? 

3.  Can  one  use  the  Thiele  algorithm  to  construct  or 
evaluate  this  rational  function? 

4.  can  one  tel^  a  priori,  when  we  can  expect  the  Thiele 
algorithm,  the  e-algorithm,  or  the  Pade  method  to  produce  an 
accurate  low  degree  rational  approximation? 

5.  Does  the  error  of  this  rational  function  compare 
favorably  with  the  error  of  the  best  possible  rational  approx¬ 
imation  of  the  same  degree? 

Although  we  cannot  give  an  affirmative  answer  to  the 
above  questions  in  all  cases,  we  shall  describe  classes  of 
analytic  functions  which  house  nearly  all  of  the  cases  en¬ 
countered  by  the  author  in  applications .  and  for  which  the 
answer  to  each  of  the  above  questions  is  "Yes". 

We  shall  develop  a  class  of  rational  approximations  for 


interpolation  over  [-1,1],  [0,oo]  and  [-oo,»].  These  rational 

approximations  share  many  of  the  features  of  the  sine  methods 
summarized  in  [17].  The  interpolation  points  of  these  ra- 
tionals  are  the  same  as  the  sine  points  and  the  classes  of 
functions  which  low  degree  rationals  approximate  accurately 
are  the  same  as  the  classes  which  the  sine  functions  approxi¬ 
mate  accurately.  Indeed  the  error  bounds  for  e.g.  approximation 
on  [-1,1]  of  functions  analytic  on  the  unit  disc  are  the 
same  as  the  sine  bounds,  i. e.  rationals  have  the  same  opti¬ 
mality  properties  as  sine  methods.  In  using  rationals  instead 
of  sine  functions,  we  lose  many  of  the  simple  relations  that 
sine  functions  satisfy,  such  as  orthogonality  and  ease  of 
getting  other  formulas  such  as  quadrature,  approximation  of 
derivatives,  methods  of  solving  differential  equations,  etc. 

However,  the  well-known  rational  function  algorithms  of 

Thiele  [21]  (the  0  -algorithm)  ,  Pade  [  12  ],  and  Shanks  [15],  Wynn  [23]  (the 

e-algorithm  )  all  share  simple  methods  of  prediction,  which 

the  sine  functions  do  not  appear  to  possess.  This  paper 

provides  an  understanding  in  that  it  enables  us  to  tell  a- 

priori,  when  we  can  expect  these  algorithms  to  work  effectively. 

The  spaces  of  functions  for  which  the  rationals  provide 
accurate  approximations  are  described  precisely  in  Sec.  2  of 
the  paper.  One  such  space,  consists,  roughtly  of  functions 
analytic  on  an  interval  with  possible  singularities  at  end- 


points  of  an  interval,  such  that  the  functions  are  of  class 
Lip^fa  >  0)  on  the  interval,  i. e. ,  functions  which  one  encoun¬ 
ters  in  nearly  all  cases  in  applications. 

The  rational  approximations  of  this  paper  have  the 
following  additional  features. 

(a)  There  are  no  poles  on  the  interval  of  approximation. 

(b)  The  rational  functions  are  linear  in  f,  the  func¬ 
tion  that  is  being  approximated. 

(c)  They  are  nearly  optimal.  More  precisely,  we  prove 
the  following  result: 

THEOREM  1.1:  Let  1  <  p  <  «,  let  p*  =  p/(p-l),  let  U  denote 

the  unit  disc  in  the  complex  plane,  let  g  be  _in  the  Hardy 

2 

space  H  (U) ,  and  let  f(z)  =  (1  -  z  )q(z).  Let  F  denote  the 

p  n 

space  of  polynomials  of  degree  n  and  set 


d.  1)  6n  =  inf 

g,eP2N+2  ’  ^^N+l  geVU) 

Then  there  exist  positive  constants  C^,C 

only  on  p  such  that  for  all  N  >  NQ, 

(1.2)  C1N"1exp{-TT(|r)1/2}  £  SNi  C2N 


sup  sup  | f (x) — 

>  !|g’!p=i  -1<X<1 

2  an<i  N0  depending 


1/(2P')expt-n(Jr)1/2) 


u 

0 


The  rational  functions  of  this  paper  are  of  the  form  y,/c 
in  (1.1)  and  they  approximate  f  on  [-1,1]  to  within  an  error 
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bounded  by  the  right  hand  side  of  (1.2). 

A  typical  approximation  result  of  the  present  paper  is 
the  following: 


THEOREM  1.2:  Let  f  and  g  satis  fy  the  conditions  of  Thm. 
1.1,  and  define  z^  and  B(z)  by 


(1.3) 


z-z  . 
_ 2. 


8,11  =  a'!!|  "i-"  >-Y  ' 


If  q  _is  selected  by  the  formula 


(1.4) 


q  =  exp(-TT(-^)  L//2} 


Then,  for  all  integers  N  >  0, 

(1.5) 


f(z  )B(z) 

max  I  f  (x)  -  E  .  „  - - rtT7 - T 

-l^l  ^-N  (z-2j)B  (Zj) 


.  „l/(2p'  )  ,  ,_L.l/2, 

<i  C2N  exp[-n(^7)  j. 


where  depends  only  on  p. 

Due  to  their  simplicity  of  construction  and  approximation 
properties,  the  rational  function  approximations  of  this  paper 
play  ^  $i milftr  role  as  the  interpolation  polynomials  obtained 
by  interpolation  at  the  zeros  of  the  Chebyshev  polynomials 
play  for  polynomial  approximation.  In  order  to  describe  this 


role  effectively,  we  return  first  to  the  case  of  Fourier  series 
Let  R  >  1,  and  let  A  denote  the  annular  region  in  the 


complex  plane  <E,  AR  =  [w  e  (E:  R  <  |w|  <  RJ,  let  F  be 

analytic  in  A  ,  and  let  C.  be  determined  from 
R  j 

/i  ^  1  2N  .  l9k  ~i;j9k.  _  2krr 

(1.6)  c^.  -  2N+1  zk=Q  F(e  e  ;  2N+1' 


Then 


(1.7) 


0^2tt 


|F(e10)  -  £j=_N  c.elj9|  =  0(R-N). 


The  bound  on  the  right  hand  side  is  essentially  best  possible 

with  regards  to  order,  in  that  the  number  R  cannot  be 

replaced  by  a  larger  positive  number  regardless  of  how  the 

c.  are  chosen. 

1 

In  (1.7)  we  now  consider  only  those  functions  F  for 
which  F (z)  «  F ( 1/z ) .  Then  we  obtain  a  cosine  polynomial  appro 
ximation.  The  mapping 


(1.8)  z  =  |(w  +  ”) 


transforms  the  annulus  A  onto  the  ellipse  E  with  foci  at 

R  R 

z  =  +1  and  sum  of  semi-axes  equal  to  R.  conversely,  if 
f (z)  is  analytic  and  uniformly  bounded  in  E  then  we  can  use 
(1.8)  to  get  a  new  function  F(w)  analytic  in  AR  with  Fourier 

series  expansion  _ _  c^e^9,  and  where  c_k  =  c^.  If 

Tn(x)  =  cos  N9,  where  x  =  cos  6,  and  x^  =  cos ( (2k-l)-^]  ,  then 
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N  f(WX)  -N 

(1.9)  max  | f  (x)  -  T^Tyr 77  )  I  =  °<R  >> 

-l^i  k  1  (x  V  W 

-N 

where  once  again,  the  R  in  the  0(R  )  bound  on  the  right 

hand  side  cannot  be  replaced  by  a  larger  number,  regardless 
of  how  a  polynomial  of  degree  N-l  is  chosen  to  approximate 
f  on  [-1,1]. 


Hence,  instead  of  finding  the  polynomial  which  best 
approximates  f  on  [-1,1],  it  is  much  easier  to  use  the 
Chebyshev  polynomial  for  which  the  interpolation  points  x^ 
are  known  explicitly  to  get  an  approximation  which  is  nearly 
as  good.  The  rational  functions  of  this  paper  share  this 
feature. 


Notice  that  for  the  case  of  polynomial  approximation 
above,  we  required  a  knowledge  of  a  region  of  analytic ity  of 
f,  a  property  which  we  can  usually  determine  a  priori  in 
applications.  Once  we  have  identified  such  an  ellipse  E 

R 

(resp.  an  annulus  A  )  we  cam  be  certain  that  polynomial  (rsp. 
Fourier  polynomial)  approximation  will  work  very  well  on 
[-1,1]  (rsp.  on  [0,2tt]).  From  the  point  of  view  of  approxi¬ 
mation  in  applications,  we  can  thus  identify  functions 
analytic  in  E  (rsp.  A  )  with  polynomials  (rsp.  Fourier 
polynomials)  since  they  can  be  very  accurately  approximated 
with  polynomials  (rsp.  Fourier  polynomials)  of  low  degree. 

Unfortunately  there  is  a  drastic  change  in  the  rate  of 
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convergence  of  polynomial  approximation  in  the  case  when 
the  function  to  be  approximated  has  a  singularity  on  the 
interval  of  approximation,  a  situation  often  encountered  in 
applications.  For  example,  if  0  <  a  <  1,  we  have 


(1.  10)  max  |  (l-x2)a  -  p  (x)  | 

-l£x£l  N 


where  p„(x)  is  any  polynomial  of  degree  N  in  x  and  C 
N 

is  a  constant  independent  of  N.  If  a  =  1/4  we  would  have 
6 

to  take  N  >  10  to  get  3  places  of  accuracy. 

While  for  practical  purposes  functions  with  singularities 

on  the  interval  of  approximation  cannot  be  identified  with 

polynomials,  there  is,  nevertheless,  a  class  of  functions 

with  singularities  on  the  interval  of  approximation  which  we 

describe  in  this  paper,  and  which  lends  itself  to  accurate 

rational  approximation.  Such  a  class  includes  the  functions 

which  we  can  accurately  approximate  with  polynomials  and 

for  practical  purposes,  we  can  identify  this  class  with 

rational  functions.  For  example,  by  Theorem  1.2  above, 

tAefe  exists  <*■ 

given  an  integer  N  >  0, rational  function  p_„  (x)/q_„  . (x) 

/  2N+2  2N+1 

with  p„  ,  of  degree  N+l  in  x  and  q„  of  degree  N  in  x, 

N+l  N 

such  that 


.  .  ,,  2  .a  P2N+2(X) 

(1.11)  max  |  (1-x  )  -  - - 7-t 

"2N+1  x 


i  x  „  Ma/2  ,  acN.  1/2 . 

I  <  C  N  exp(~tr(— )  ]  . 


We  remark  that  by  identifying  classes  of  functions 
which  can  be  approximated  accurately  by  rational  functions, 

we  are  identifying  classes  of  functions  for  which  we  can 

^  ' 
expect  the  Pade  method,  the  Thiele  algorithm  or  the  Epi/- 

lon  algorithm  work  well.  We  shall  later  in  this  paper 

illustrate  this.  For  example,  we  would  be  able  to  tell  a 

priori,  that  the  Pade  method  used  in  [2]  may  be  expected  to 

be  accurate. 

Another  practically  important  use  of  rational  functions 

is  in  analytic  continuation.  For  sake  of  illustration,  let 

us  momentarily  return  to  the  class  of  functions  analytic 

and  bounded  in  the  ellipse  E  described  above.  Let  us  assume 

R 

that  f  is  known  on  [-1,1],  and  that  we  want  to  evaluate  f 

at  the  point  ■“  +  ^{R  +  -^)  in  the  ellipse.  This  can  be  done 

by  means  of  the  polynomial  in  (1.9),  the  rate  of  convergence 

N 

of  the  error  zero  being  0(p  ),  where 

p  =  [  (a+l)/2  +  va+(a2-3)/4]/R,  a  =  (R  +  R_1)/2.  On  the 
other  hand,  if  j*N(x)  is  any  polynomial  approximation  to 

2  eg  11, 

f(x)  =  c  +  ( 1— x  )  on  and  we  want  to  approximate 

f ( 1 )  =  c  by  evaluating  PN(1),  then  we  may  expect  [ f ( 1 ) —  p^ ( 1 ) ] 
to  converge  to  zero  very  slowly,  indeed,  too  slowly  to  be 
of  any  practical  use.  Since  however,  we  may  identify 

2  d 

f(x)  =  c  +  (1-x  )  with  a  rational  function  for  practical 
purposes,  we  can  accurately  evaluate  f(l)  via  a  rational 


function,  by  using  values  of  x  on  only. 

As  a  more  sophisticated  example,  let  u  =  u(x,y)  be 
harmonic  in  the  right-half  plane,  and  assume  that  u(0+,y) 
is  of  class  Lip  (a  >0)  on  a  neighborhood  of  y  =  0.  It 
follows,  then,  that  u(x,0)  is  analytic  and  bounded  on  a 
sector  with  vertex  at  the  origin  and  of  class  Lip^  on  [0,A] 
where  A  >  1  is  arbitrary.  That  is,  for  practical  purposes, 
we  may  identify  u(x,0)  with  a  rational  function,  and  we 
can  accurately  approximate  u(0,0)  via  a  low  degree  rational 
function,  by  using  values  of  u(x,0)  on  (e.g.)  the  interval 
[1,A]. 

In  the  cases  when  the  condition  of  accurate  approxima¬ 
tion  are  satisfied,  it  is  thus  possible  to  do  accurate  analy¬ 
tic  continuation  all  the  way  to  the  boundary  of  analyticity, 
via  a  relatively  low  degree  rational  function. 

The  Lip  property  of  the  function  to  be  approximated 
ct 

is  important  from  the  point  of  view  of  applications;  if  f 
approaches  zero  too  slowly  in  a  neighborhood  of  a  singularity, 
then  it  is  necessary  to  choose  the  degree  of  the  rational 
function  to  be  very  large,  in  order  to  achieve  a  desired 
accuracy.  For  example,  for  rational  approximation  on  [0,1], 
if  f(x)  -  f(l)  ~  c/[log(l-x)  ]  as  x  1  then  it  is  just 

as  difficult  to  approximate  f  on  [0,1]  by  a  rational 

2  ot 

funv-wx  n  ;  it  is  to  approximate  (1-x  )  on  [-1,1]  by  a 


polynomial.  We  remark  however,  that  this  difficulty  can 

often  be  remedied  by  means  of  a  transformation.  For  example, 

if  we  set  x  =  l-exp(-z  ^)  ,  we  get  f(l-exp(-z  1))  -  f(l) 
a 

~  cz  ,  z  -*•  0,  and  we  can  now  approximate  the  new  function  of 
z  defined  as  the  interval  [0,°o]  by  a  rational  function. 

We  mention  that  a  rational  function  was  previously 
constructed  by  the  author  [18]  of  the  same  degree  as  that  in 
(1.5)  for  approximating  f  on  [-1,1],  and  moreover  the  error 
bound  in  [18]  is  the  same  as  that  on  the  right  hand  side 
of  (1.5).  However,  whereas  the  interpolation  points  in  [18] 
are  the  points 

(1.12)  =  k1/2  sn[(2j-l)K/(2N)  yk] 


the  evaluation  of  the  * .  is  more  difficult  than  the  evaluation 

a3 


of  the  z j  in  (1.3). 


The  same  points  z_.  defined  in  (1.3)  were  also  used  by 


Peaceman  and  Rachford  [13]  to  approximate  the  points  in 


(1.12)  in  their  alternating  direction  iterative  method  for 
obtaining  approximate  solutions  to  elliptic  partial  differ¬ 
ential  equations. 

For  many  problems  of  rational  approximation  one  does 


not  have  analyticity  in  the  unit  disc  U,  but  rather  in  a 


smaller  region  7^  (see  Figure  2.1  in  Sec.  2)  and  we  have 


therefore  also  considered  this  case.  Although  cur  error 


bounds  in  this  case  are  not  as  small  as  the  sine  bounds,  we 
suspect  that  the  rationals  of  this  paper  are  as  accurate  as 
the  sine  approximations  for  the  same  problem,  and  that  it 
may  be  possible  to  improve  the  bounds  of  this  paper  in  the 
case  when  0  <  d  <  rr/2. 

Notice  that  if  N  is  replaced  by  4N  in  the  rational 
function  in  (1.5"  )  then  every  second  interpolation  point  in 
the  "4N"  rational  case  is  the  same  as  the  interpolation  point 
in  the  "N"  rational  case,  and  the  "4N"  rational  approximation 
has  roughly  twice  as  many  correct  significant  figures  as  the 
"N"  rational  approximation.  This  result  is  of  practical 
value,  particularly  when  a  user  is  unable  to  determine  a 
region  (Figure  2.1). 

Let  us  now  briefly  describe  the  layout  of  the  paper. 

In  Sec.  2  we  give  precise  statements  and  proofs  of  the 
results  (a),  (b) ,  (c)  and  (d)  stated  at  the  beginning  of  this 
section.  These  proofs  would  ordinarily  be  lengthy,  and  for 
this  reason  some  of  the  details  are  carried  out  in  appendix 
A  and  B. 

In  Sec.  3  we  illustrate  connections  of  the  results  with 
the  well-known  approximation  algorithms,  the  Thiele,  or 
p-algorithm,  the  Epsilon  algorithm  (c -algorithm)  and  the 
Pade  method.  In  view  of  the  results  of  Sec.  2,  we  are  able 
to  determine  a  priori,  when  we  can  expect  these  algorithms 


to  work. 


In  Sec.  4,  we  prove  Theorem  1.1  above.  While  the  exact 
optimal  rate  of  convergence  of  rational  approximation  is  not 
known,  we  conjecture  that,  in  the  notation  of  (1.1), 

(1.13)  in-  sup  sup  |f(x)  -  |  ^ 

5  aePN  ?eHp(u)  •  !!gllp=1  -1<x<1 

as  N  -*>  a>. 

In  Appendix  A  we  study  bounds  on  rational  functions 
related  to  (l.S).  The  Jacobi  theta  functions  turn  out  to  be 
most  convenient  for  this  purpose,  since,  while  it  is  possible 
to  obtain  similar  results  via  the  approximate  integration 
of  the  function  F(z;t)  =  t  ^  log  j  (z+t)/ (z- 1)  |  over  0  £  t  £  °o 
via  the  trapezoidal  (resp.  midordinate)  rule  evaluated  at  the 
points  (rsp.  ) ,  j  =  0,+l,+2,...,  and  using  the 

concavity  of  this  function  for  (fixed  z*(6,°°))  as  a  function 
of  t,  it  is  possible  to  achieve  more  accurate  error  bounds 
via  the  theta  functions  since  it  is  possible  to  get  exact 
bounds  via  known  properties  of  the  theta  functions.  However, 
while  we  use  elliptic  functions  to  obtain  cur  results,  the 
final  results  are  independent  of  elliptic  functions. 

In  Appendix  B  we  obtain  accurate  bounds  on  contour 
integrals  encountered  in  the  proofs  in  Sec.  2. 

We  close  this  section  with  a  few  historical  remarks. 
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Stietjes  [20  ]  seems  to  be  the  first  to  identify  classes  of 
function  which  may  be  approximated  via  truncated  forms  of 
continued  fractions.  These  functions  are  expressible  in 
the  form 


and  the  continued  fraction  expression  obtainable  via  this 
representation  converges  uniformly  in  any  closed  region  of 
the  complex  plane  which  does  not  contain  the  interval  [a,b] 
(see  [7]).  Unfortunately  given  a  function  F  it  is  not 
possible  to  easily  check  in  applications  whether  or  not  F 
has  a  representation  of  the  form  (1.14). 

In  [6]  Gautschi  gives  an  excellent  summary  of  the  use 
of  rational  functions  in  numerical  analysis.  It  has  long 
been  suspected  and  verified  in  ad  hoc  cases  that  a  rational 
function  can  do  a  better  job  of  approximation  in  a  neighbor¬ 
hood  of  a  singularity  than  a  polynomial.  The  first  quantita¬ 
tive  result  as  to  exactly  how  much  better  a  rational  function 
can  be  than  a  polynomial  was  obtained  by  Newman  [10,11]. 
Renewed  interest  has  developed  in  rational  functions  since 
Newman's  result.  The  error  of  the  rational  functions  of  this 


paper  have  the  0(e 


-cn 


1/2 


rate  of  convergence  which  was 


obtained  first  by  Newman.  Saff  and  Vargr.  (see  [14,22])  have 


obtained  many  beautiful  results  demonstrating  the  superior 
power  of  rational  functions.  Also  of  interest  is  the  idea  of 
Ganelius  [5] ,  for  using  the  Greens  function  of  a  region  of 
analyticity  to  obtain  rational  approximations;  indeed  the 
rational  functions  of  this  paper  have  this  property.  For 
the  case  of  rational  approximation  on  a  finite  or  semi¬ 
infinite  interval,  the  poles  of  the  rational  functions  of 
this  paper  lie  on  the  real  line  outside  of  the  interval  as  is 
the  case  for  best  approximation  of  Stietjes  transforms — see 


Borwein  [3] . 
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2.  RATIONAL  APPROXIMATION  WITH  ERROR  BOUNDS 


This  section  contains  the  main  approximation  theorem 
of  the  paper.  While  the  proofs  are  complete,  we  use  results 
derived  in  the  appendix  in  order  to  shorten  these  proofs. 

At  the  outset  we  cover  in  detail,  the  case  of  rational 
approximation  on  the  interval  [-1,1].  These  results  are 
then  extended  to  the  case  of  rational  approximation  with  a 
rational  function  of  a  variable  g  over  a  contour  P, 
where  ^  is  a  one-one  transformation  of  P  onto  the  interval 
[-1,1].  We  then  use  this  generalized  result  to  obtain  two 
types  of  rational  approximations  over  the  interval  [0,a>]  and 
one  over  [-°°,°°]. 

Rational  approximation  on  the  interval  [-1,1]  is  effect¬ 
ive  when  applied  to  a  certain  class  of  functions  that  is 
analytic  in  the  region 

(2.1)  e  C:  |ar^ir5l  <  0  <  d  ^  tt, 

where  E  denotes  the  complex  plane.  In  the  case  when 

d  =  tt/2  ,  is  the  unit  disc.  When  0  <  d  <  tt/2,  is  the 
d  a 

intersection  of  the  two  discs 

(2.2)  2  =  fc  €  E:  jc  +  i  cot  d|  <;  |csc  dj} 

,  is  the  union  of  the  two 
d 


while  if  tt/2  ^  d  <  tt 
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discs  (2.2).  We  shall  obtain  rational  approximations  for 
functions  analytic  in  3)^,  and  also  for  functions  analytic  on 
the  regions. 

(2.3)  $d  =  (z  e  p:  | arg  z|  <  d], 

3 

(2.4)  $d  =  [z  €  E:  jarg  sinh  z|  <  d), 

(2.5)  $d  =  [z  e  E:  | Zm  z|  <  d). 

These  regions  are  illustrated  in  Figs.  2. 1-2.4. 

2.1  Rational  Approximation  on  [-1,1], 

We  describe  two  typical  situations  of  rational  approxi¬ 
mation  on  the  interval  [-1,1].  The  conditions  in  the  first 

case  are  of  theoretical  interest,  particularly  when  $.  is 

d 

the  unit  disc,  while  the  conditions  in  the  second  case  can 
be  readily  tested  in  applications. 

Given  f  analytic  in  3)^,  we  define  F  by 

(2.6)  F(C)  = 

i-c 

Assumption  2.1a?  Let  f  be  analytic  in  35^,  let  F  be 
defined  by  (2.6),  and  for  some  p  in  the  range  1  <  p  <  oo,  let 

(2.7)  !|F|!  .  f  ,|F(C)|P|<a:|)1/p'  < 

0 


Let  a  be  either  0  or  let  0  <  q  <  1,  let  be  defined  by 


<2-8>  Cj  - 

3  q3  CT+1 


and  corresponding  to  a  positive  integer  N,  let  B(f)  be 
defined  by 

c-c< 


N 


(2.9)  B(C)  =  n"=_ 


j=-N+2cr  1-CjC 


Theorem  2.  la«  Let  Assumption  2. la  be  satisfied,  let 


0  <  d  tt/2,  and  define  T)  by 


(1-5  )B(§)  f  (£  .  ) 

(2.10)  1|(S)  =  f (?)  -  Ej=.N+2CT  - 1- 


(5-Cj)  (l-C-)B' (Cj) 


where  B  ( Q )  _is  defined  in  (2.9).  _If  q  _is  selected  by  the 
formula 


(2.11)  q  =  expf-rr^)1^2}. 


where  p'  =  p/(p»i) ,  then  there  exists  a  constant  C  depend in« 
only  on  p,  such  that  for  -1  <,  %  <>  1. 


(2.12)  |71(f )  |  ^  CN1/(2p,)exp{-d(|r)1/2},!F!!p. 

Proof;  It  is  readily  seen  that  71(5")  in  (2.10)  also  has 
the  representation 


T|  (?) 


2rri 


_ f(C)dc 

(c-?)  u-c2)b(c) 


(2. 13) 


Now,  by  Thm.  A. 10  (i),  the  function  ( B ( ^ ) | ” 1  satisfies  the 
inequality 

,  tt  CH-d)  +  e  (3,  q) 

(2.14)  max  | B { C )  |  ^  expf - ; - } 


,  1 
log  - 

q 


where,  by  (A.  36),  we  have,  for  0  <  d  tt/2  , 


(2.15)  ,  exp(--^)  ^  _n_ 

log  —  log  —  log  —  2ed' 

q  q  q 


Hence,  setting 


(2.16)  A  *  exP^2d  e 


taking  absolute  values  of  each  term  in  (2.13)  and  using  (2.6) 
we  get 


(2.17) 

(2. 18) 


|T)(;)|  <;  0^2Hb<OIr 


.  1 
log  — 

q 


} 


a®.. 


um 

lc-si 


|dC 


2  "(“-d)  t' 

£  Ml-S  )  |B(g)  lexp(-- -■■-1-}[G(pl,d,g)]P  || F | 


log  ~ 

q 


where  (2.18)  is  obtained  after  applying  Holder's  inequality 
to  the  integral  on  the  right  hand  side  of  (2.17)  and  then 
using  the  notation  of  Eq.  (B.18)  of  Appendix  B.  By 
applying  Thm.  B. 1  to  bound  G(p',d,|),  we  therefore  get 


(2.19) 


1,(5)  |  i  A  21/P’(sin  a,  <p'-l>/(2p' )  [-LiEU-L—j  1/p' 

_  /£_  \  /P  \ 


(2 . 19  cont.  ) 


2  l/o'  "(2d) 

•  <l-5^)1/p  |B(;)|exp{ — — H. 

log  - 

q 

We  shall  now  attempt  to  obtain  a  uniform  bound  on 
(1-§^)^P  | B  ( f )  | .  To  this  end,  let  us  select  a  number 
such  that 

(2.20)  0  <  qN+1/2"CT  ?1  <  1 

and  let  us  define  I(§^)  by 

S,*1  l-S, 

(2.21)  KljlMx,— ). 

Then,  by  Thm.  (A. 10) 

(2.22)  max  |B  (§)  |  <;  b(N,q,  %  ) 

feKSj^) 

where 

2 

(2.23)  b(N,q,§1)  =  2  exp{-  - 33 - -  (l-qN/§  )  } . 

2  log  - 

q 

Since  |B(f)|  ^  1  on  [-1,1],  we  have 

if  5  e  K?L) 
if  56[-l,l]-l(g1) 


2  1/p'  i  ^ 

(2.24)  (l-z  )  /P  |B(?)  |  ^  <  4. 

■  ( - - — )  ^P ' 

-  v  2 

U+5,) 


We  now  define  q  by  (2.11),  5.  by 


to  get  (2.12),  m  which  we  may  take  C  to  be 


(2.26)  C  =  2Vp'  (sin  d)  2p'  f~  )P'  expfn [  (2ed) ~  K  (2p‘  )  2]  ) 

t(p/z)  t(f/2) 

Remark  ( i )  :  If  d  =  rr/2,  the  result  (2.12)  remains  the  same, 
except  that  the  constant  C  may  be  taken  to  be 


(2.27)  C  =  21+1/p'exp{TT(2p*  )"1/2}f- 


Remark  (ii)  :  From  (2.12),  the  dominant  rate  of  convergence 
of  the  error  to  zero  as  N  ■>  »  is 

(2.28)  rL  =  exp[-d  (~“) 1/2 } . 

The  corresponding  rate  for  sine  approximation  is 

(2.29)  r  =  exp{-(^N)1/2}. 

^  £r 

This  rate  r2  converges  to  zero  more  rapidly  than  r^  if 
0  <  d  <  tt/2,  and  r^  =  r^  if  d  =  tt/2.  The  author  originally 
expected  these  rates  to  be  the  same  for  all  d  in  the  range 
0  <  d  <,  rr/2,  but  we  have  so  far  been  unable  to  obtain  the 
sharper  bound  of  (2.29). 

Remark  (iii) :  It  would  be  interesting  to  extend  Thm.  2.1a 
to  the  case  rr/2  <  d  <,  n*  Thm.  A.  10  applies  to  this  case; 
it  shows  moreover  that  we  must  more  carefully  bound  the 
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contour  integral  (2.13),  taking  advantage  of  the  fact  that 
whereas  |B(£)|  ^  is  close  to  1  in  neighborhoods  of  Q  =  +1. 
it  is  close  to  2  expf 1/q^  °n  Part  3®^  outs*-de 

of  the  neighborhoods  of  +1.  We  have  not  done  this,  due  to 
the  added  complexity  of  the  proofs.  However,  such  an 
approach  may  also  lead  to  an  improvement  of  the  bound  (2.12) 


for  the  case  when  0  <  d  <  tt/2  . 


Remark  (iv) ;  The  Assumption  2.1a  is  not  satisfied  if  f(*) 

does  not  approach  zero  as  |  +1.  However,  if  g  is 

analytic  and  bounded  in  3)^  and  is  of  class  Lip^fS^), 

where  3D,  denotes  the  closure  of  3),,  then  f  defined  by 

d  d 

(2.30)  f(E)  =g(?)  -^g(-D-^g(l) 

satisfies  Assumption  2.1a,  for  p  <  l/(l-a).  After  obtairtng 
a  rational  approximation  for  f  we  add  the  linear  term 
( 1— 5 ) g  (—  1 )  +  ^(l+Cjgfl)  to  get  a  rational  approximation 
for  g. 

Rather  than  testing  whether  or  not  Assumption  2.1a  is 
satisfied,  it  is  often  simpler  in  practice  to  check  whether 
or  not  the  following  assumption  is  satisfied.  All  of  the 
above  remarks,  albeit  after  some  obvious  modifications,  apply 
also  to  this  case. 


Assumption  2.1b;  Let  f  be  analytic  in  3),  and  let 


(2.31)  |f (C> |  ^  c | (1+c)  (1-C )  | 


I 

for  all  r  _in  where  c,  a  and  3  are  positive  numbers , 
and  0  <  a  <  1,  0  <  3  <  1. 

Theorem  2, lb :  Let  Assumption  2 . lb  satisfied ,  and  let 

0  <  d  ^  tt/2.  Let  y  =  min(a,f3),  5  =  max(a,<3),  and  correspond¬ 
ing  to  a  positive  integer  n,  let  q  bg  defined  by 

(2.32)  q  =  expC-T-r^Yn)-1^2} . 

If  y  =  a,  let  M  and  N  be  defined  by 

M  »  n 

(2.33) 

N  =  [fnl 

while  if  y  =  let  M  and  N  be  defined  by 
M  =  [|n] 

(2.34) 

N  =  n. 


Let  q ^  be  defined  by  (2.8),  and  let  B  ( q)  and  T)  (§)  be  defined  by 


BO  =  nj*-M+2o-  7^ 


(2.35) 


jst7 
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Lfl 


t* 


(2.36) 


N  (1-x  )B(5)f (;  ) 

n(?)  =  f(s)  -  ^j=:-M+2CT  ~ — 77  2 - 1- 


(?-Cj)(l-Cj)  B'(r.) 


Then  there  exists  _a  constant  C  depending  only  on  c,  a, 

3  and  d,  such  that  for  all  *  _in  [-1,1]  , 

(2.37)  |T)(5)|  £  C  n5/2  exp{-d  (2yn) 1/2  } . 

Proof :  The  identity  (2.13)  now  also  defines  the  T)  (f)  of 
(2.36),  provided  that  B(£)  is  defined  by  (2.35).  Taking 
absolute  values,  replacing  |f(£) |  by  the  right  hand  side  of 
(2.31),  and  using  (2. 14)- (2. 16) ,  we  get 

2  TT(?~d) 

(2.38)  |71(?)|  £  Ac(l-?  )jB(g)|H(a,3,d,g)expC . 

log  - 

where  H(a,6,d,p)  is  defined  by  Eq.  (B. 19)  of  Appendix  B. 
Using  the  result  of  Thm.  B.2,  we  now  get 

g  ft  n(?-d) 

(2.39)  |Tl(5)  |  £  A  C  C^l+f)  (l-g)P|B(c)  |exp(  ■  '  ’  1  ) , 

log  - 

where  the  constant  depends  only  on  a,  3  and  d. 

We  now  proceed  to  obtain  a  uniform  bound  on 

a  s 

(1+?)  (I-?)  |B(§)[.  To  this  end  let  and  =  ^  be  selected 

such  that 


(2.40) 


0  <  qM+1/2-a  £  5i  <  j. 


o  <  qN+l/2^  ^  ;2  <  1 


and  define  I(£^+£2) 


V1 


(2.41)  I(?L,c2)  =  Cff*  ■— j  £  g£ 


Then,  by  Thm.  A. 10, 


(2.42)  max  ) B ( % ) |  £  b (M, N, §  , §  ) 

?€i(c1,52) 

where 


(2.43)  b(M,N,Px,§2)  ^2exp{- 


M  N 


2  log  - 

q 


1[1  “  2?1  "  2?2]K 


We  therefore  have,  since  |B(x) |  £  1  if  -1  £  x  £  1, 


20+0-3 


‘r1 


^  -q+i 


(2.44)  (l+;)a(l-5)®|B(5)  I  £  {  b(M,H,?1(|2)  lfseK^.S^I 


l  ,4+20  0 

N  -2 


Now  q  is  given  by  (2.32)  and  we  consider  first  the  case 
when  y  =  a  =  min(a,0).  In  this  case  M  =  n,  N  =  [an/3].  If 
we  furthermore  select  and  §  by  the  formulas  5^  =  qMM^^ , 
'2  =  q^N1^ ,  then  we  get  the  uniform  bound  (2.37). 

The  argument  in  the  case  when  y  =  2  is  similar,  and  we 


omit  it. 


2.2  Rational  Approximation  on  a  Contour 


Let  S  be  a  simply  connected  domain  with  boundary 
935,  let  a  and  b  (b  ^  a)  belong  to  as,  and  let  9  be 
a  conformal  map  of  S  onto  the  region  S^  (Eq.  (2.1))  such 
that  cp(a)  =  -1,  cp  (b )  =  1.  Let  $  =  p  ^  denote  the  inverse 
map,  and  define  r  by 


(2.45)  T  -C*(C> s  -1  £  C  4  !)• 


Let  2^  denote  the  point 


(2.46)  z  =  *(Cj), 


where  the  £  are  defined  as  in  (2.8).  Let  .€  c.e  analytic 
in  S  and  define  y  by 


(2.47)  y(z)  = 


l-cp(z)' 


Let  58  (z)  be  defined  by 

N  C5  (z  ) 

(2.48)  5B(z)  =  rL_  M4-  -7— - rfr 

3  =-M+2cr  l~CjCp(z) 

where  a  is  defined  as  in  (2.8). 

Let  one  of  the  following  two  assumptions  be  satisfied: 


Assumption  2.2a:  Let  5  be  defined  by  (2.47),  and  for  some 


I 
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(2.49)  j|y||  s  lim  infH^  f  |y  (z )  | P '  (z )  |  | dz  j  ) 1//p  < 

"  UC 


Assumption  2.2b:  Let  f  be  analytic  in  3),  and  for  all 
z  e  3),  let 

(2.50)  |f(z)  j  ^  Ajl-cp(z)  |3|  1+cp  (z )  \P 


where  A,  a  and  3  are  numbers  such  that  A  >  0,  0  <  a  <  1, 
0  <  3  <  1. 


Proceeding  as  in  the  previous  section,  we  set 


(2.51)  Tl(x)  = 


f  (z)fa1  (z)dz 


€  r. 


a®  [cp  (z)-cp  (z)  ]  [l-cp(z)  J  8  (z ) 


to  get 

[  l~cp  (x) 2  ]  a  (x)cp  1  (z  ) 

(2.52)  T)  (x)  =  f(x)  -  mj~ - ^ - 1—  • 

:,='M+2a  b{x)-cj][i-c-]a,(2j) 

The  sum  on  the  right  hand  side  is  a  rational  function  in 
the  variable  m(x). 


Theorem  2.2a:  Let  Assumption  2.2a  be  satisfied,  let  0  <  d  £  rr/2 , 
let  M  =  N  _in  (2.48)  and  (2.52),  and  let  r  be  defined  by 
(2.47).  Lf  q  _is  defined  by 

(2.53)  q  =  expf-Tr(^“) 1//2} 


where  p '  =  p/ ( p- 1 ) ,  then  there  exists  a  constant  C 


depending  only  on  p  such  that  for  all  x  _on  r 

(2.54)  |T](x)  j  £  CN1/2p'exp(-d(^7)1/2)!j^'!  . 

1  p  P 

Proof;  If  we  set  <q(z)  =  q  in  (2.51),  we  get  (2.7)  with 
y(’Jf(2))  =  F  (z )  ;  if  we  set  «j(f)  =  x  and  0(z)  =  j  in  (2.51) 
we  get  (2.13).  Hence  the  proof  is  identical  to  the  proof  of 
Thm.  2.1a. 

The  proof  of  the  following  theorem  is  also  similar  to 
the  proof  of  Thm.  2.1b,  in  view  of  the  above  remarks. 

Theorem  2.2b;  Let  Assumption  2.2b  be  satisfied,  and  let 


o  <  d  ^  rr/2.  Let  y  =  min  (a,  S)  .  5  =  max  (a,  3)  ,  and  correspond 
ing  to  a  positive  integer  n,  let  q  be  defined  by 


(2.55) 

q  =  exp(— it (2yn)  7  }. 

If  y  =  a. 

let  M  and  N  be  defined  by 

M  =  n 

(2.56) 

N  =  [fn] 

while  if  - 

/  =  6,  let  M  and  N  be  defined  by 
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Let  £  be  defined  by  (2.8),  z  ^  b^  (2.46),  F  b^  (2.45)  and 
3J(z)  and  T)(tf) 

N  sp  (z  ) 

(2.58)  .(,)  =  nj,-M,2g 


[l-cp(x)  ]8(x)f (Zj)cp'  (z  ) 

(2.59)  T|(x)  =  f (x)  -  E?  - - -  . 

3-M+2°  [#(x)-cj)a-Cj)»,(*j) 

Then  there  exists  a  constant  C  depending  only  on  a  and  B 
such  that  for  all  x  e  r. 


(2.60)  |  T)  (x)  |  ^  C  n6/2  expf-d  (2yn) 1/2  ) . 


2.3  Rational  Approximation  on  [0,oo];  the  Non-oscillatory  Case. 


In  this  case,  the  region  s  of  analyticity  is  the 
2 

sector  of  Eq.  (2.3).  The  function  qj  and  the  inverse 
function  of  the  previous  section  are 

(2.61)  z  =  *(c)  =  ^  «  C  =  cp(z)  = 

Hence  corresponding  to  the  points  rv  in  (2.8),  the  points  z_. 


(2.62)  z.  =  q- 


(c  =  0  or  1/2). 


Thus  the  product  (2.48)  becomes 


b(z)  =  z-q~!~g 

3  M+2<y  z+q3-a 


(2.63) 


Corresponding  to  f  analytic  in  the  function  7  in  (2.47) 
is  now  defined  by 

(2.64)  ?(z)  =  t^-ffz); 

the  condition  (2.49)  thus  becomes 

(2.65)  lim  (~  I  .  |*(2)  |P  <  cot  1  <  P  <  °>> 

t-0+2"^(t)  |1«|2 

2 

where  ®  (t)  =  [z  €  E:  |arg(2-T)|  <  d],  while  the  condition 
(2.50)  becomes 

(2.66)  jf(z) |  ^  A|l+2|‘a“0|z|P. 


The  error  T],  i.e.  ,  the  difference  between  f  and  the 
rational  function  thus  becomes 


(2.67) 


T](x) 


f(x) 


«B(x) 
1  +  x 


-M+2p 


(l+qj~^f  (z.) 
q^Vz.l?1  (z..  ) 


2 

It  thus  follows  that  if  f  is  analytic  in  S)d>  and 

(2.65)  is  satisfied  we  take  M  =  N  in  (2.63)  and  (2.67), 

and  choose  q  as  in  (2.53)  to  bound  T|  on  [0,<»]  according 

2 

to  (2.54).  If  f  is  analytic  in  2D,  and  (2.66)  is  satisfied, 

a 

then  by  choosing  q  as  in  (2.55),  we  bound  T|  on  [0,oo] 


according  to  (2.60). 
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2.4  Rational  Approximation  on  [0,°o]  ;  the  Oscillatory  Case. 


In  this  case  the  region  3)  of  analyticity  is  the  region 
5)^  of  Eq.  (2.4).  The  mapping  functions  $  and  ijr  are 
defined  by 


1  +  s  inh  ^  (■—■£ ) 


.  sinh  z  -  1  .  .  xl-r 

(2.68)  o(z)  -  sinh  z  +  i>  *<C)  ~  .  -l  i+£ 

1  -  sinh  (~ — *) 

i-c 


We  see  therefore  from  (2.63)  that  gj(z)  takes  the  form 


(2.69)  S8 


.  .  N  sinh  z  -  q3  CT 

<z)  *  ^  j=-M+2(r  ~  7h 

J  sinh  z  +  q 


i.e.,  the  points  z..  are  given  by 


(2.70)  Zj  =  sinh  1(q^  ff)  =  logfqJ  u  +  J l+qZJ  CT). 


If  f  is  analytic  in  33^,  the  corresponding  function  y 
is  defined  by 


<2-7i»  *<z>  =i¥rfi¥I- f(z)- 


The  condition  (2.49)  becomes 


(2.72)  ilT||  -  f  3  1 5 (z)  j P  |2— °3h  Sl2  lds|)P  <  «, 

^  t-*0+  n  J  (t)  |l+sinh  z| 

1  <  P  <  °°> 


(t)  =  [z  €  E:  | arg  sinh(z-T)j  <  d] , 


where 


and  the  following  condition  is  equivalent  to  (2.5 0) 


(2.73)  |  f  (z)  j  <,  A 


f  (z)  I  ^  Alz  j3|e  32  I ,  z  e  3^ 


The  difference.  T|>  between  f  and  the  rational  function 
becomes 


(2.74)  T)(x)  =  f(x) 


j-o  f  (z  .  )/  1+q2  (j-c) 


.-jinh_x  n  i±slLI  -  --i::-:-- - 

l+sinh  j=-M+2c  j-c r  (sinh  x-z.)8'(z.) 

q  j  3 


Either  Thm.  2.2a  or  Thm  2.2b  may  now  be  stated  to 
bound  T](x)  on  [0,<x>]. 


2.5  Rational  Approximation  on  [-00,00]  Via  a  Rational  Function 

^  x 
of  e  . 


The  region  3)  of  analyticity  is  the  strip  35^  defined 
in  Eq.  (2.5).  The  function  m  and  $  become 


(2.75)  ®(2)  = 


e2-l 

eZ+l 


•i(C)  *  logr1^* 


The  product  8(z)  now  takes  the  form 


(2.76)  s(z)  =  n 


N  e2-q^  g 

j=-M+2c  z .  j-q’ 
e  +q 


i.e.  the  points  z_.  are  given  by 


(2.77)  z.  =  ( j-<r)  log  q. 


Let  f  be  analytic  in  3)^.  Then  the  corresponding  function 
y  of  (2.47)  is  defined  by 

z  2 

(2.78)  y(z)  =  1  V  f (z). 

4e2 


The  condition  (2.49)  now  becomes 


(2.79) 


=  lim  (  ; 

T-»d  w335(t) 


(1+eV 


I  dz  I ) 1//p  £  oo,  p 


The  following  condition  is  equivalent  to  (2.50): 


(2.80) 


f(z) 


fA  exp  (a  Re  z)  , 
A  exp (-(3  Re  z). 


4 

z  e  Re  z  <  0 

d 

4 

z  e  3).,  Re  z  >  0. 
a 


The  difference  between  f  and  the  rational  function  now 
becomes 


(2.81) 


Tl(x) 


f  (x) 


x 


1+e 


-M+2a 


(W  <7)f(zi) 
(eX-^  Isb'  (z  j  ) 


We  may  now  bound  T|  ( z )  on  (-«,<»)  via  either  Thm.  2.2a  or  Thm. 
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3.  IMPLICATIONS  AND  APPLICATIONS. 

In  this  section  we  study  the  connection  of  the  results 

of  the  previous  section  with  the  Thiele  algorithm,  the 

/ 

epsilon  algorithm,  and  the  Pade  method. 

3.1  The  Thiele  Algorithm. 

The  Thiele,  or  p  algorithm  for  interpolating  f 


m  4-  1  distinct 

points  xQ 

i  X«  >  •  •  •  >  x 
x  m 

is  described  as  follows 

Define  by 

r  'l  ’ 

f(x.) 

j  =  0,1, 

.  .  .  ,TV> 

i)  <j  = 

i 

X.  -X. 

■  It1 _ 1 

j+1  3 
p  0  -p  0 

j  =  0,1, 

.  .  .  , m- 1 

j  _ 

X  -X. 

X+l  X 

.  j  +  1 

j  =  0 , 1 , . .  .  ,  m-  c 

^  pi  " 

j+1  j 

+  p  i-2 

i  “  2,3,..., m. 

pi-rpi-i 


Then  the  rational  function  r{x)  which  interpolates  the  data 
fXj , f (Xj ) } is  given  by  the  continued  fraction  representation 


(3.2) 


r(x) 


+ 


x  -  x.  i 

X  I 

— .  i 

0  0 

p2  “  p0 


+  • 


.  + 


The  function  r(x)  has  the  form 


if  m  =  2n,  where  pn  and  are  polynomials  of  degree  n  m 
x,  and  it  has  the  form 


(3.4) 


Vi (x) 


r(x)  - 


q*(x) 


if  m  =  2n  +  1.  Furthermore,  if  m  =  2n,  then 


(3.5) 


r  (x)  = 


On  n-1 

p_  x  +C.X  +...+  C 

K2n _ 1 _ n 

n  n-1 

x  +  d,x  +. . , +  d 
1  n 


That  is 


(3 •  6)  P2n  =  lim  r  ( x) , 
X^+oo 


so  that  the  Thiele  algorithm  provides  an  excellent  method  of 
carrying  out  analytic  continuation. 

For  example,  if  f  is  analytic  and  bounded  in  the 


region  5)  of  Eq.  (2.3)  and  if  f  is  of  class  Lip  (a  >  0) 
u  a 

[Xq,<»],  where  xQ  J>  0,  then  we  may  effectively  use  the 
Thiele  algorithm  to  evaluate  f(«)  via  the  use  of  a  few 
values  of  f(x)?  for  finite  x.  Indeed,  this  has  been  done 
recently  in  an  ultrasonic  tomography  algorithm  (19 ] . 


3.2.  Evaluation  of  the  Rationals  of  the  previous  Section 


via  the  Thiele  Algorithm. 


Let  B  denote  the  family  of  polynomials  of  degree  £  n 


and  consider  the  evaluation  of  the  rational  function 


for  p 

6  E  ,  q  e  P  ,  where 

a  =  0  or  1, 

and  such  that 

n"  ^n+c;  n+cr 

(3.8) 

fr<*2k>  -  f'x2k)’ 

k  =  0,1,... 

,n 

< 

l  r(x2k-l’  =  ”  ’ 

k  =  1,2,... 

,  n+cr . 

Then 

(3.9)  pOO-T^j 


can  be  evaluated  via  the  Thiele  algorithm,  using  the 

=  l/r(x^),  k  =  0,1,...,  2n  +  1  +  cr- 
Eq.  (3.9)  then  yields  r(x)  =  1/p (x) .  In  general,  there  seems 
to  be  no  guarantee  that  the  p  algorithm  will  always  work; 
however,  interlacing  the  zero  and  non-zero  values  of  p 
in  the  above  fashion  has  worked,  in  our  experience.  Since 
all  of  the  poles  of  r(x)  have  been  pre-determined,  there 
are  no  unwanted  poles. 

Let  us  next  consider  the  evaluation  of  the  rational 
function  in  (2.10)  for  the  case  when  cr  =  0,  i.e.  , 


2n  +  1  +  a  values  p (x^) 


r  (x) 


(1-x  )B(x)f(;_.) 
(x-Cj) (l-Cj)B’ (Cj> 


(3. 10) 


=  4=^ . 

a1  ±1 


b (x)  =  n 


N  X~€i 

-N  1-r  .x 
WD 


Since  rQ  =  0,  B(x)  =  p^+1  (x)/q2N(x)  ,  where 

P2N+1  £  E2N+1’  q2«  5  E2N-  HSnCe  r,X)  *  P2N+2(x,/q2N(lt) 
2 

and  it  has  (1-x  )  as  a  factor.  Hence 


(3.  11) 


p  (x)  = 


1-x 

r(x) 


P™(*> 

q2N<X>' 


P2N,q2N  €  P2N  } 


is  completely  determined  by  the  4N  +  1  values 

N+k 


(X2k,p(X2k))  (C-N+kJ 


'si 


f {^-N+k) 


), 


(X2k-l"(X2k-l^  = 


(1/<:-N+k-l’0)’ 


(i/c.N+k.°). 


k  =  0, 1, . . . , 2N 


k  1 , 2  ,  .  .  .  ,  K 

k  =  N+l, ...  ,2 


and  may  thus  be  evaluated  via  the  Thiele  algorithm,  as  above. 
The  rational  function  r(x)  may  then  be  computed  via  (3.11), 
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On  the  other  hand,  if  S ^  is  defined  by 


(3.16)  Sj  =  c.t,1 


n 


then  e2k  yields  a  Pade  approximation  [23]  i.e.. 


(3. 17) 


n  pm+k- 1 ( 11  * 
®2k  =  qk(H) 


The  results  of  Sec.  2  of  this  paper  together  with  the 
representations  (3.13)  or  (3.16)  tell  us  when  we  may  expect 
the  approximations  (3.15)  or  (3.17)  to  be  accurate,  when 
applied  to  a  function  f. 

For  example,  if  f(x)  -  L  (L  =  f  (00)  )  satisfies  any  of 
the  conditions  in  sections  2.4  or  2.5,  then  by  taking 
si  =  f(3k)>  we  see  from  (3.13)  and  (3.15)  that  e°  will 
converge  rapidly  to  L  =  f(<»). 

The  representation  (3. 16)  may  be  considered  to  be  an 
interpolation  of  an  analytic  function  f(z),  where 

<3-18)  w-viov* 

This  representation  then  shows  that  f(z)  -  L  (L  =  f(»))  may 
be  assumed  to  satisfy  any  of  the  conditions  of  Secs.  2.4  or 
2.5,  provided^  -e . 

(3.19)  |Ckxk|  =  0(e‘3k),  a  >  0. 

Hence  in  the  case  when  (3.19)  is  satisfied,  we  may  expect 
(3.17)  to  converge  rapidly  to  f(<»)  =  z^_Qckxk  as  n  - 


A  LOWER  BOUND  ON  THE  ERROR  OF  BEST  APPROXIMATION. 


Let  u  denote  the  unit  disc  m  the  complex  plane,  i.e 


3)  in  the  notation  of  Eq.  (2.1).  That  is 
tt/2 


(4.1)  U  =  IQ  £  <C:  \q  j  <  1}. 


Let  f  satisfy  the  conditions  of  Thm.  2.1a.  Then  F 

defined  by  (2.6)  is  in  H  (U) ,  that  is,  F  is  analytic  in 

P 

U,  and 


(4.2) 


=  lim 
r->l 


"-‘ir 


|F(rei8)|pd9)1/p 


Let  S  denote  the  family  of  all  such  functions  F  such 
that  ||F||^  £  1.  Then,  by  Anderssen  [1]  we  have 

(4.3)  C.  expf-rr  (^7) 1//2)  =  inf  sup|  [K  ( 5 )  - — 2—  ■  . 

P2N'*2N+1  PSS  *  -l  92H+1' 

where  p'  =  p/(p-l)  and  where  is  a  positive  constant  de¬ 
pending  only  on  p,  and  p^N  and  ^n+I  denote  polynomials  of 
degree  2N  and  2N  +  1  respectively. 

Now  if  T](p)  is  defined  by  (2.10)  then  the  terms  in 
square  brackets  on  the  right  hand  side  of  (4.3)  is  just 
T)  (f  )/(l-f2) .  Hence 

C.  exp  { — tt  (~7 ) 1/2  }  <;  if  d?|. 


(4.4) 


1 
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We  now  split  the  integral  on  the  right  hand  side  into 
an  integral  over  (-6*5),  5  >  0  plus  an  integral  over 
[-1,1]  -  (-6,6). 

Then 


(4.5) 


■^^2  d?|  £  max  1 7]  (f)p  (l-x~)  ^dx 
J-s  1-1  5€[-i,l?  v  -  6 


-  l09(-rr-J) 


while  from  (2.fJ)  we  have 

(4.6)  ||"  JLlSl 

^[-1, 1]- (-6,6)  1-5 


dD  £  C. 


f\l  -,V1/P 


6 


^  c2  r  (i-D'1/pdi 


=  C2P? (1-6) 


P' 


where  c2  is  a  constant  depending  on  on  p. 

Hence,  by  (4.4),  (4.5)  and  (4.6),  we  get 


H 


(4.7) 


lInlL  > 


expf-n(|r)1/2)  -  2p’U-6)l/P. 


2  ^ 

Taking  1  -  6  =  e  ,  and  combining  with  Thm.  2.1a  we  find 

that  there  exists  an  NQ  ^  0  such  that  if  N  >  N  ,  then 

there  are  constants  C,  and  C.  such  that 

3  4 


(4.8) 


<; 


inf  sup  sup  |f(?) 

P2N+2,q2N+l  FSS  "1<?<1 


cy/^P'^expf-.^)1/2)  . 


q2N+l(^ 


These  inequalities  show  that  while  the  exact  lower  bound 
on  the  left  hand  side  of  (4.8)  is  not  known,  the  results  of 
this  paper  are  in  the  right  "ballpark"  with  respect  to  their 
accuracy. 

We  therefore  conclude  with  a  problem:  Given  f 

2 

analytic  in  U,  F  €  Hp(U),  where  F(z)  =  f(z)/(l-z  ), 

®  * 

and  given  N,  can  a  rational  approx  i  t*aho*  of 

f  on  [-1,1]  which  is  linear  in  F  be  as  accurate  as  the 
best  rational  approximation  to  f  of  the  form  Here 

p,  and  q T  are  polynomials  of  degree  at  most  N. 

*  XT  XT 
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Appendix  A.  BIASCHKE  PRODUCT  ESTIMATES. 


We  shall  see  in  what  follows  that  the  four  rational 
approximation  results  referred  to  at  the  outset  of  Sec.  1 
are  simply  related  via  elementary  conformal  transformations. 
For  example,  once  we  have  a  rational  function  approximation 
in  w  over  [-1,1]  we  can  readily  obtain  one  in  z  over 
(0,a>)  via  the  transformation 


(A.  1) 


w 


z-1 

z+1 


(«  z 


1+w 

1-w 


which  is  a  conformal  map  of  the  right  half  plane  onto  the 

unit  disc,  and  which  transforms  p^+^ (w)/qn<w)  into  a  rational 

function  P  ,.(z)/Q  _(z).  Similarly,  starting  with  a 

n+l  n+1 

rational  function  of  z  for  approximation  over  [0,«>]  we 
can  use  the  second  equation  in  (A. 1)  to  get  a  rational 
approximation  in  w  over  [-1,1]. 

It  is  most  convenient  from  the  point  of  view  of  using 
known  results  for  Jacobi  theta  functions  to  first  consider 
approximation  over  [0,»].  For  purposes  of  interpolation 
at  the  points  ,  j  =  -N,...,N,  0  <  q  <  1  it  is  natural  to 
start  with  the  product 


$*  (z)  =  Z~<L. 

V  -D-N  z+qD' 


0  <  Z  <  00. 


Unfortunately  this  product  does  not  have  a  limit  as  N  »,  since 
the  product  changes  sign  with  N  as  N  increases.  However 


the  alternate  form 
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(A.  Z) 


,  _  Zz±  _N  1_-  gUz+l/z)  +  ,g2j 

%(Z)  “  z+1  nj  =  l  .  j  ,  .  .  .  2  j 

J  1  +  q  (z+l/z)  +  q 


has  the  same  zeros  and  poles  as  $*,  and  moreover,  the  product 


converges  as  N  ■*  =°.  Let  us  set 
(A.  3) 

and  let  us  define  m(q)  by 


slz  a)  _  Zzl  r*  g3  (z+i/z_L..+_.q2  j 

z+1  3-1  x  +  q3 (z+l/z)  +  q23 


lA.^)  m(q)  =  sup  |$(z,q) 

0<Z<oo 


Next,  let  us  use  the  standard  notation  for  elliptic 
functions,  for  0  <  k  <  1, 

rV. 


(A.£1 


/-  u  =  u  (k)  = 


dt 


J 


«  w  =  sn(u;k) 


0  I  ,1  U2\  ,1  !  2  2i 

J  (1-t  )  (1-k  t  ) 


cn(u;k)  =  J 1-sn2 (u?k) ,  -K  <  u  <  K 


K  =  K(k)  =  u  ( 1 ) 


=  J  1-k2 


K'  =  Kfk^ 


q  =  expf-Ti  K'/K] 


qx  =  expf-rr  K/K'] 


We  then  prove 

Lemma  A. 1 :  If  m (q)  is  defined  by  (A. 5 ) .  then 


,  2, 
m(q  ) 


(A.  4) 


Proof:  Let  us  introduce  the  Jacobi  theta  functions 


Ss(u)  =  (■ 1//6sin  v  n”  ,  (l-2q2^cosiv+q4^)v 
S  J~L 


(A.  7) 


9  (u)  =  ( - r-¥^6cos  v  n°°  ,  (l+2q2:*cos  2v+q4^  )\ 

k  D-l 


Then,  by  (A. 3)  and  (A. 7)  we  have 


(A. 8) 


2iv  2  1/2  Vu> 

5(e  -q  !  =  lki  Tm- 

c 


Now,  using  Eqs.  (16.36.3)  and  (16.3.3)  of  [9]  we  have 


(A.  9) 


»<e2iv,q2)  =  ik^2 

^  1  cn(u;k) 


However,  we  are  interested  in  bounding  $(z,q)  on  the 
interval  0  £  z  £  <».  This  means  that  we  want  to  find  the 
maximum  value  of  |  $  (e2i  ^_iV^  ,q)  | ,  for  -  oo  <;  v  However 

by  Eqs.  (16.20.1)  and  (16.20.2)  of  [9]  we  have 


(A.  10) 


sn  (-iu;k) 
cn (-iu?k) 


-isn  (u?^) 


But  for  u  real 


(A. 11) 


so  that 


max  |sn(u;k^) 

-a^U^oo 


=  1, 


l+qJ  '-(z+l/z)  +  q  J 

for  2N  -  point  rational  approximation  at  the  points 
<r(j~i/2),  j  =  1,2, ...  ,N. 

Letting  N  <».  we  set 

(A.  14)  T(z,q)  .  n"  1  ~  ^  ^  <z+l/z)  +  g2-i~1 


(A. 15)  n(q)  = 


sup  |y(z,q) 
0<Z<°° 


We  then  prove 


Lemma  /{^_2 :  _If  n  (q)  _is  defined  by  (A.  15)  ,  then 


(A.  16)  n(q2)  =  k  Y2- 


Proof:  in  this  case  we  use  the  theta  functions 


.2  2 

.  ,  1,1/12  oo 

!a(u)  =  (ir;>  ni: 


oo  2  i  —  1  2 

+  2q  cos  2v  +  q 


(A.  17) 


Vu>  *  <■ 


Js _ ,1/12  «  ,,  ,  2 j-1 

4 '  -  2q  cos  2v  +  q 


and  proceed  similarly  as  the  proof  of  Lemma  2.1,  using 
equations  in  [9], 


Lemma  A.  3  :  Let  k^  and  q  _be  defined  as  in  Eqs ,  (4  .g) . 

Then 


(4.18)  kj-/2  ^  2  expf-  — “ - -}. 

4  log  - 

q 

Proof;  The  following  series  expansion  taken  from  [8,  p.  378] 
converges  for  0  <;  k^  £  1: 

(A.19)  q1/4  =  ) 1/2  [  1  +  2(|)2  +  15  (“ ) 4  +  150  (~)  6  +...]. 

Since  all  of  the  terms  of  the  series  in  square  brackets  are 
positive,  (4.19)  taken  together  with  (iK.S)  implies  that 


(A.  20)  \  £  2  q^  . 


Since,  however,  the  last  two  equations  in  (A. S')  yield 


(A.  21)  (log  ^)(log  ~ )=  n2; 


we  find,  after  substituting  this  result  into  (A. 20)  and 
using  the  fact  that  0  <  q  <  1,  0  <  q  <  1,  we  obtain  (A. 18) 


Lemma  A.  4 :  The  functions  $  of  (A.  4)  and  y  _of  (A.  14 )  satisfN 


(A.  22)  sup  |  $  (z,q)  I  2  expf- 


(A. 23) 


sup  |y(z,q)j  £  2  exp{-  - r). 

0<z<co  2  log  — 

q 

Proof:  As  a  consequence  of  Lemmas  A. 1,  A. 2  and  A. 3,  both 
2  2 

functions  $(z,q  )  and  Y(z,q  )  are  bounded  on  [0,»]  by 

2  1  2 
2  expf-rr  /(4  log  — )].  Hence,  replacing  q  by  q  we  get 

(A. 22)  and  (A. 23). 


We  remark  that  while  we  used  elliptic  functions  to 
derive  (^.22)  and  (W.23)  these  inequalities  do  not  involve 
elliptic  functions. 


Lemma  A.  5 :  Let  a  be  either  0  _or  1/2. 
N+l/2-a 


2  I  ^  ^  q 


(A. 24) 


,  and  if 


z+q~ 


N  j  =W+1  z_qj-. 


then 


2  N+l/2-a 

(A.  25)  1 1  |  £  exp{ - n - -  H - - - ). 

4  log  —  ® 

q 


(b)  If  0  <  |z|  ^  q'N  1/2+ff,  and  if 


(A. 26) 


-1  j-a 

z _ +  q 


M  “j=M+l  -1  j  — or 

z  -  qJ 


then 


/a  -37 \ 


i  -  i  / 


— n -  „  rt+i/2-q. 


(a)  If 


Proof:  The  proof  of  (A. 27)  is  almost  identical  to  the  proof 
of  (A. 25),  and  hence  we  shall  only  prove  (A. 25).  We  have 


(A. 28) 

But 

(A. 29) 

Hence 

(A.  30) 

This  is 

Lemma  A 
Then 

(A. 31) 


iog|tN|  =  Re  z”=0  ^  ^K+1<4f>2n+1 

”  2  a>  cp  f.2n+l 

*  En=0  2n+l  E j  =N+1  \z\  } 

.  oo  2  oo  ,q^  c ,  2n+l 

^  En=0  2n+l  Ej=N+l  ?  * 

.  N+l/2-(j  „  ,  n+1/2 

_  r°°  2  ■«.  _ _ 2n+l  _g _ 1_ 

_ r\  •*>  _  .  i  K  —  / 


"11=0  2n+l  5 


l-q 


2n+l* 


_ _ 1  ^  ■  1 

l-q2n+1  2  sinh[ (n+“) log  ~]  (2n+l)log 


_  S 


N+1/2- 


n+1/2-  <T 

5 _ 1 _  »  _ L_ 

§  log  ~  n=0  (2n+l) 2  5  log 


just  the  logarithmic  form  of  (A. 25). 


_6:  Let  |z|  >0,  — rr  <  9  <  tt>  and  let  z  -  Izle 


ie 


I  t  -  TT(f  -  (9  |)  - 

vJ  0 


Proof:  Splitting  the  integral  as  an  integral  from  0  to 

z  plus  an  integral  from  iz|  to  oo  we  have 


•°  H  L  a  im 


(A.  32 )  Re  jQ  -  log  —  dt  =  *.j  ^ 


_  2  -(2n+l)iS  °°  2 

=  Re  - 2e  =  n - rcos  (2n+l)  9 

n— u  .  « % 2.  ri— 0  .  _  _ .  2 

(2n+l)  (2n+l) 


Similarly, 

CO 

(A. 33)  Re 

I  - 

Hence 

00 

■1 

(A.  34)  \ 

J  n 

~  log 

=°  (2n+l)  2 


cos  (2n+l)  e. 


=0  2 
( 2  n+ 1 ) 


cos  (2n+l)  a 


=  "<?  -  |9|) 


from  Fourier  series. 

Corollary  A,  7 :  z  satisfies  the  conditions  of  Lemma 

/\ .  &  and  if  9  -  (t:  t  =  pe1®,  cp  fixed.  0  p  <£  <»} ,  then 

(A- 35)  jffli  £ 

Proof :  The  proof  is  similar  to  the  proof  of  Lemma  A. 5,  by 

proceeding  along  the  lines  of  equations  (A.  3/)  t<?  (A.  3  4.), 

dropping  the  real  parts,  and  replacing  t  by  p  and 
z  by  |  z  | . 


Lemma  A,  8 :  Let  0  <  J  9  |  <  rr ,  let  z  =  j  z  |  e~  ,  and  define 
e (9,q)  by 
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r 


(A. 36)  e (B,q)  = 


2 

_  H_ 
2 


exp )  exp (-2-n(TT-  |9D  } 
log  -  log  - 

- 2_  +  - L<I - 


l-exp(-^4)  1-exp  (-^(TT-l6i)) 

-  1 _ _  ■»  -L 


log  — 

q 


log  — 

q 


If  <7  _is  either  0  or  1/2,  then 


(A. 37) 


1  00 


W<2  -  |  9  | )  -  log  -f  z-  log  I  I  <;  «<e,q> 


j+CT 


q  3 


z-q 


Proof ;  Without  loss  of  generality,  let  us  fix  9  in  the 
range  0  <  9  <  tt.  The  function 


(A. 38)  f(t)  =7  log  ?-+-  ^ 

z  -  tq‘ 


is  then  an  analytic  function  of  t  in  the  sector 

9  -  tt  <  arg  t  <  0,  it  is  absolutely  integrable  along  any 

ray  t  =  1 1 1  e  ^ ,  0^  |tj  9  -  tt  <[  tp  <[  9,  and  moreover 

for  t  in  this  sector  of  analyticity.  f(t)  =  0(1)  as  t  -*  0, 
-2 

f(t)  =  0(t  )  as  t  -m»,  Hence  if  we  can  apply  Thm.  4.3 

of  [17],  Lemma  A. 6  and  Corollary  A. 7  to  get  the  result  r^.3 6) 
for  0  <  9  <  tt.  The  proof  for  the  case  when  -rr  <  9  <  0 
is  similar. 

i  9 

Corollary  A.  9 ;  If  z  =  j z  |  e  ,  where  0<  |  9  |  <rr,  if_  e  ( 9 , q ) 
is  defined  by  (A. 36 )  and  if  g  is  either  0  or  1/2 ,  then 


52 


(A.  39)  exp(- 


tt(£-|  q|  )-e(e,q) 


CO  ,  Z-K 


1  ^  nj=“OT  I  j+C 

log  —  z-q 

q 


i  exp{- 


TT  (-?-  I  9  I  )+€  (9,q) 


Theorem  A. 10;  Let  z  =  jz|e  ,  let  €<o,q)  be  defined  by 
(A. 36)  and  let  be  either  0  or  1/2. 

(i)  if  0  <  |9|  then* 

m  „  ,  _3 —  cr  tt (~~  |  o  |  )  +  e(3,q) 

(A-40)  <  -<■ — I^J — 1 


(ii)  if  tt/2  <  9 1  <  tt,  let  §  and  r  be  selected  such 


th3t 


(A.  41)  qN+1/2-c  ^  ?  <  c  ^  q-(Wl/2-a) 


if  |  z  |  lies  is  the  ranee  ?  <,  |  z  |  ^  c  >  thej? 

2  N-K^-<r  M4 

j -cr  n(f-|9|)+e  o..q)+^-^q  /?+q 

<A-«>  "  eXP' - - 

(iii)  If  g  =  tt.  and  z  lies  in  the  range  %  <  \z\  <,  r, 
where  %  and  Q  satisfy  (2.41),  then 


*If  1 9  I  =  the  left-hand  side  of  (A. 40)  is  identically  1. 


w  Im 


(A-43)  3j=M+2ff  ■  j-c '  *  2  eXP'- 

z-qJ 


log  — 

q 


or  equivalently ,  if  ^  <,  z  <,  £,  then 


N  z-o^-CT  ”  2  4  ^ 

(A'44)  S=-M+2a!-^l  *  2  - ; - 1 

z+q  log  ^ 

o 


Proof:  (i)  We  have 


2  2  N+~-<t  M+~-cr 

2  /?  +  q  2  c 


(A  45)  nN  \S±£ll  I  -  r,ra  ,z± 

nj=-M+2a'z_qj-cl  Uj=-~'2_( 


|2±a_,  ,  »  i^±a_|n  is=ai 

1  z_q^  1  “3  — 1  z_q3-a  1  J  j<-M+2a,  j>N '  j 


However  since  |  G  |  =  |  arg  z\  <>  rr/2  each  term  (z-q3  CT)/(z+q ^ 
of  the  product  on  the  right  hand  side  of  (4. 45)  is  at  most 
1  in  magnitude.  Hence  (A. 40)  then  follows  from  (A. 39). 

(ii)  This  is  also  a  consequence  of  (A. 45),  (A. 39), 

and  Lemma  A. 5. 

(iii)  Both  |$(z,q)j  (by  taking  a  =  0)  and  |y(z,q)| 

(by  taking  c  =  1/2)  are  expressible  by  means  of  the  infinite 

product  nj__oo|(z+q^  ~)/(z-q^  a)  j  .  Hence  the  inequalities 


(A. 43)  and  (A. 44)  follow  from  Lemmas  A. 4  and  A. 5. 
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f  <6,  p.  50] 
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(B.  8 ) 

(B.  9) 

(B.  10) 

(B.  11) 

Some 

(B.13) 

(B. 14) 

(B. 15) 


^ 

F(a,b;c;z)  =  (1-z)  F (c-a ; c-b ; c; z ) 


[8,  P-  47] 


F  (a,b;c;z)  =  ,  ~~lj!  F  (a, b ;a+b-c+l ;  1-z ) 

T (c-a) r (c-b) 

+  ( 1— z )  C  3  b  (c-a ,  c-b ;  c-a-b+1;  1-z ) 


18,  P-  47] 


F(a,b;c;z)  =  <l-z>‘*  ^^f^P(a,c-bta-b+l,^j) 

+  (l-z)_b  ,L(C)  ^ L§rJ?.) p (k  c-a*b-a+l- — — ) 

(  J  r(a)r(c-b)F(  '  a,b  a+1'l-z) 

\b>  p.  48] 


F(a,b;c;z)  =  (1-z)  aF (a, c-b;cr~-) 

1-z 


=  (l-z)"bF(b,c-a;c?-~) 

[8,  p.  47] 


gamma  function  identies  will  also  be  convenient: 

r(2)r(z  +  |)  =  n  1/221'2zr(2z) 


r(z)rd'Z) 


_ 31 _ 

sin  ttz 


+  z)r(-2  -  z) 


_ rc _ 

cos  (rrz) ' 


Let  be  defined  by  t 


The  Integrals  to  be  Estimated 


We  shall  estimate  the  integral 


(B.16)  G(p’ ,d,x)  *  \  lC“*l  P  ldC|  1  <  P*  <  °°> 


-1  <  x  <  1, 


0  <  d  £  tt/2 


for  l<p'  <«,  -1  <  x  <  1,  and  0  <  d  ^  tt/2,  and  we  shall 
estimate  the  integral 

(B.  17)  H(a,S,d,x)  =~  f  )l+c,ct-1,l-c,3-1,c-x|-1!dC| 


for  0  <  a  <  1,  0  <  3  <  1,  -l<x<  1,  and  0  <  d  <;  rr/2. 


Theorem  B.  1 :  Let  G(p' ,d,x)  be  defined  by  (B. 16) .  Then 
there  exists  a  constant  C  depending  only  on  p*  and  d 


such  that 
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over  parts  of  two  circles.  If  we  replace  the  integrals 
over  these  parts  by  the  integrals  over  the  whole  circles, 
we  find  that 


(B.23) 

(B.a*) 

(B  .Z5) 

(Compare  (B. H  ) ,  (B.  )).  The  .  inequality  (B./S) 

now  follows  by  proceeding  as  in  (B.<2<2)  -  (B.22.). 

Let  H(a,8,d,x)  be  defined  by  (B.j*7j,  and  let  us  set 


G(p’,d,x)  <,  i  j  C“v  x2+cot2d  |  P  |d;| 

TTw  j  £  I  =csc  d 

=  2 (esc  d)^  P  ~  1  I £-sin  d J x2+cot2d j  P 

»  2  (esc  d)^  ^  G(sin  dl/x^+cot^d)  - 


(B.2W  H!x)=^j  1 1+C  I9'1  IC-^r1  ldC  |  • 

Icl-l 


2  i 

Then,  setting  q  -  e  ,  0  <;  g  tt,  we  get 

j 

3-1 

cos  gsm 


1  rr/2 

(B.2L7)  H (x)  =— -r—7  f  cos3~18sin3"13fH-— ~X  ,sin23]  ~1/2de 

nU"X)  (1-x)2 


(B.2S) 


(B  .13) 


_  2a*3~2  F(2^  r(2)  1  8  a+g  4x 

"(1-x)  _,3+S.  F(2’2?  2  7‘  2 ’ 


r(*f) 


(1-X)' 


by  (B.j).  Then  (B.$)  yields 


,  .  23+d~2  r(  2  )r(2)  .1  8  1+x  2 

H(X)  *  n(l-x)  r  .a»9-l)  F  2 ' 2  7  2  ;  (l-x}  ] 


,3+8-2 


r  (— )  r  (— - ) 

...  ,3-1..  .-3  2  _. 3+8-1  3  3+1  .1+x. 

(1+x)  (1-x)  - J F  ( — - — ?  (T~x) 


r(^) 


On  the  other  hand,  replacing  sin  Q  by  cos  ;p  in  (B.^7) 
and  then  using  (B .4),  we  get 


H  (x)  = 


1o+b-2  r(f)r(^)r[i  g  3.6  ^)2) 


tt (1+x)  i  (q+S-l)  2 ’ 2  2  1+x' 


(B.50) 


.a+@-2 


r  (— )  r  ^ ) 

— (i-x^-Ni.x)-9 l,pa=iifp, 

r  (2 ) 


Let  us  next  transform  (B./7)  by  means  of  the  transfor¬ 
mations 


,  w-1  u— 1 

<B-3" 


to  get 


CB.  3-d)  H(x) 


o+0-l  f08  a-l,,  2,  2~^,  2.  2  ,  -1/2 

=  - - v  M  w  (1+w  )  (w  +u  )  c 

2tt  J0 


,a-l,  ,6, 

r(-r-)rCr) 


2a+0-3  *  x  2  ' 1  '2 ' _ ,1  6  3-a  2. 

~  (1+U)  "r':oT6-r~F(2-2?~7U ) 

n  2 


,o+S-3 


O  1-0 

..  ..  a-l1"  2  r  2  _,a+S-l  o  a+1  2, 

— (1+M)U  — X - F(— —'2  — !U  > 


rtf  > 


from  (B.2.9J*  and 
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(B.33) 


.  .  2a+3~3  1+u  r(2)F(  2  ]  1  a  3-3  -2 

H  X  tt  u  r(a±6=l)  F  2  *  2  *  2  ,U  * 


2a+S"3  -0-1  r(2)r(  2  l  a+6-1  3  3+1  -2 

+  - u  (l+K) - ; - F  ( — ~ — ), 

"  r(^) 


from  (B.3c). 


The  representation  (B.2£)  is  convenient  on  the  interval 

-1  <  x  ^  0,  while  the  representation  (B.^3)  is  convenient 

on  the  interval  0  <;  x  <  1. 

1+x  2 

Now  if  -1  <  x  £  0,  (r—-)  <,  1;  by  inspection  of  the 

power  series  of  the  hypergeometric  function  in  (B.Z^)  and 

recalling  that  0  <  a  <  1,  0  <  ^  <  1,  we  find  that  each  takes 

2  2 

on  its  i*»/ix<mum  value  if  (1  +  x)  /(I  -  x)  =  1.  Hence, 
applying  (B.£  ) ,  we  get 


r(f)r(^)r(Ir3>r(1::§;:§)  2ch3-2 

(B'3*}  H(X)^  rC^r^rf^psT" 


d-x)a-1r(f)r(1f)r(If)r(^fg)2a-f3-2 

tt2  r  ^ ) 


if  -l  <  x  £  o 


B.3S ) 


H  (x)  £ 


r(®)  2a+B-2 

r^jrc^firt^f6)  " 


(i-x)e-1r(f)r(^),^)r(^)2“+6-2 

n2r<^f> 


if  0  <|  x  <  1. 


Hence,  denoting  by  C  the  maximum  of  the  right  hand  side 
of  (B.3V)  and  (B ,.3S)  ,  we  have 


(B.4C)  H(x)  ^  C(l+x)Ct"1(l-x)3  \ 


We  now  consider  the  case  0  <  d  <  rr/2.  In  this  case, 
we  consider  the  integral  obtained  by  using  (B .3/)  in 

(B./7), 


(B.3/)  H(a,0,d,x)  = 


id 


2tt 


/ 


.  (a-l.,  ,1-a-B,  -1 .  ,  , 

w  I 1+w I  w-u  dwj 


Here,  |w 


.  at—  1 , 


1+w 


l-a-0 


is  a  decreasing  function  of  |w|  along 


the  path  of  integration,  while  |w-u|  first  increases  to 
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(u  sin  d)  and  then  decreases  again.  Hence,  since  u  sin  d 
is  the  distance  from  u  +  iO  to  the  path  of  integration,  we 


have  by  rearrangement  of  the  functions,  that 


l-a-0 


-titp  .  . 

(B .3$)  H<a,0,d,x)  <,  - U 


Jet—  1  2  2  22  2 

w  [1+w  ]  [w  +u  sin  d]  aw 

o 


=  2H(u  sin  d) 


1  +  u 
1+u  sin  d 


where  the  last  identity  was  obtained  using  (B .32.). 

Theorem  B.2:  Let  H(a,3,d,x)  be  defined  by  (B./71,  where 
0  <  a  <  1,  0  <  0  <  1,  0  <  d  tt/2.  Then  there  exists  a 
constant  C  depending  only  on  a,  3  and  d.  such  that  for 
-1  <  x  <  1, 

(B. 42 )  H(o,0,d,x)  £  C(l+x)a_1(l-x)0"1. 


•— t|  fN 
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20.  ABSTRACT  (Continue  on  rovoroo  at  do  it  nocaaaory  pad  idontify  by  block  nutnbor) 

This  paper  gives  explicit  rational  functions  for  interpolating  and  approxi¬ 
mating  functions  on  the  intervals  [-1,1]  ,  [0,~]  and  [-«,»]  .  The  rational 
functions  are  linear  in  the  function  f  to  be  approximated.  The  error  of 
approximation  of  these  rationals  is  nearly  as  small  as  the  error  of  best  rational 
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approximation  with  numerator  and  denaninator  polynomials  of  the  same  degree. 
Regions  of  analylicity  of  f  are  described,  which  make  it  possible  to  tell  a 
priori,  the  accuracy  which  we  can  expect  from  rational  approximation.  A  sample 
result  is  the  following:  Let  f  be  analytic  in  the  sector 
V  =  {z€(C:  |arg  z|  <  d  <_  tt/2}  ^  define  F  by  F(z)  =  (l  +  z)^f(z)/z  .  For  seme 
p  in  the  range  1  <  p  <  00  ,  let 


#FII  =  sup  I  ‘ 
p  1 0 1  <  d  UO 


|F(ceie)|P  — iE- 
ll  +  te1 


Let  p’  =  p/(p-l)  ,  and  corresponding  to  a  positive  integer  N  ,  let  q  be 


defined  by 


Let  B(t)  be  defined  by 


q  =  exp{-Tr[p7(2N)]i/z}  . 


B(t)  *  H  . 

j*-N  t  +  qJ 


Then  for  all  t  €  [0,°°]  , 


T 

where  C  is  a  constant  depending  only  on  p  and  d  . 
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